Phases of Wave Functions and Level Repulsion 



W.D. Heiss 

Centre for Nonlinear Studies and Department of Physics 
University of the Witwatersrand, PO Wits 2050, Johannesburg, South Africa 



On 
ON 

On 



G 



> 

m . 

(N ■ 

o ■ 
o : 

On ■ 
On ■ 

Ok 
-i— > ■ 

G ■ 

g : 
cr 



PACS: 03.65.Bz, 02.30.Dk, 84.40.-x 



Avoided level crossings are associated with exceptional 
points which are the singularities of the spectrum and eigen- 
functions, when they are considered as functions of a cou- 
pling parameter. It is shown that the wave function of one 
state changes sign but not the other, if the exceptional point 
is encircled in the complex plane. An experimental setup is 
suggested where this peculiar phase change could be observed. 



The phase of the wave function changes in a charac- 
teristic way, if a selfadjoint Hamiltonian has an energy 
degeneracy at specific values in some parameter space 
and if a loop around the point of degeneracy is described 
in parameter space [[l] . In the simple case of a real sym- 
metric Hamiltonian two parameters are needed to get a 
diabolic point. In this case Berry's phase is tt when loop- 
ing around the degeneracy in two dimensional parameter 
space. 

If a selfadjoint Hamiltonian depends only on one pa- 
rameter, its variation will in general give rise to level re- 
pulsion Q . Associated with a level repulsion is a pair of 
exceptional points || which are the points where the two 
levels actually coalesce when continued analytically into 
the complex plane of the parameter [Q. In the present 
letter we discuss the behaviour of the wave functions and 
their phases when an exceptional point is encircled. We 
suggest an experimental situation where such behaviour 
could be measured. It is distinctly different from a dia- 
bolic point where Berry's phase occurs. In particular, an 
exceptional point must not be confused with a coinciden- 
tal degeneracy of resonance states, which was considered 
in (|] a s a generalisation of Berry's phase. 

The difference between a diabolic point and an excep- 
tional point is due to the selfadjointness of the Hamilto- 
nian in the former and the lack of it in the latter case. 
Also, when continuing into the complex parameter plane 
we are faced with analytic functions of a complex vari- 
able which implies a more rigid mathematical structure. 
In the references quoted above a thorough discussion is 
given of the spectrum and eigenfunctions, when one pa- 
rameter on which the Hamiltonian depends is continued 
into the complex plane. The parameter chosen can be an 
interaction strength as considered in this paper, but other 
choices are possible ||. Generically, an A^-dimensional 
matrix problem yields N(N — 1) exceptional points. For 
an infinite dimensional problem, an infinite number can 



occur [[p, but, depending on the particular situation, they 
may accumulate at one point in the finite plane j^] ; there 
are Special cases, where no exceptional point oeeurs in 
the firfi;e plane as can be seen from the analytic solution 
for the sTfrKte^carticle motion in the>-Hulthen potential 
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FIG. 1. Level repulsion using in Eq. 
uj\ = 2, LJ2 — — 1 and (j> — n/25. 



All essential aspects of exceptional points can be illus- 
trated in a two level model on an elementary level. In 
fact, for finite or infinite dimensional problems an iso- 
lated exceptional point can be described locally by a two 
dimensional problem [ jTo| . In other words, even though 
a high or infinite dimensional problem is globally more 
complex than the two dimensional problem, we do not 
loose generality for our specific purpose when the restric- 
tion to a two dimensional problem is made. For easy 
illustration we therefore confine ourselves to the discus- 
sion of 
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This is, up to a similarity transformation, the most gen- 
eral form of a real two dimensional Hamilton matrix of 
the type Hq + XHi. The particular dependence on the 
parameter A has been chosen as it is of a nature widely 
used in physical applications. We emphasise that our 
aim is not in particular directed at a physical model that 
is describable by a two dimensional problem although 
there may exist interesting problems in our special con- 
text. The restriction has been chosen for illustration, 
while the physical application that we have in mind is an 
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infinite dimensional situation. 
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FIG. 2. Contours in the complex energy plane. The top 
illustrates contours for each level, if a closed loop is described 
in the A-plane without encircling the exceptional point. The 
bottom illustrates the energy contour produced by two (equal) 
closed loops in the A-plane, which encircle the exceptional 
point. The solid line corresponds to the first closed loop in 
the A-plane and the dotted line to the subsequent loop. The 
solid dot is the position where the contours would meet, if the 
loops would cross the exceptional point in the A-plane. 



The eigenvalues of H are given by 

ei + £2 + A(wi + wa) 



£i, 2 (A) 



±R 
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where 
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+ ( A(a;i ~ W%) f + ^A( E1 - eaXm - ^) cos 2cA . (4) 



Clearly, when <f> 
lines 



the spectrum is given by the two 



Xuji 



1,2 



which intersect at the point of degeneracy A = — (ei — 
e 2)/( w i — Wg). When the coupling between the two lev- 
els is turned on by switching on <f>, the degeneracy is 
lifted and avoided level crossing occurs as is illustrated 
in Fig.l. Now the two levels coalesce in the complex A- 
plane where R vanishes which happens at the complex 
conjugate points 
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At these points, the two levels Ef.(X) are connected by 
a square root branch point, in fact the two levels are the 
values of one analytic function on two different Riemann 
sheets. In Fig. 2a we display contours in the complex 
energy plane for each level, which are obtained if a loop is 
described in the complex A-plane, which does not encircle 
the exceptional point. Correspondingly, Fig. 2b shows the 
contours, if the exceptional point is encircled. In this case 
only a double loop in the A-plane yields a closed loop in 
the energy plane. Obviously, this connection is not of the 
type encountered at a genuine diabolic point. 

The difference has a bearing also on the scattering ma- 
trix pd[ | and on the wave functions ^i(A) and ^(A). In 
pd[ , although the notion exceptional points is not used, 
the pertinent distinction between a genuine degeneracy of 
two resonances and the analytic coalescence (exceptional 
point) of two (complex) eigenvalues is nicely discussed. 
A usual degeneracy of two resonances still gives rise to a 
simple pole in the scattering matrix or Green's function, 
while an exceptional point produces a double pole. With 
regard to the eigenfunctions, we recall that for complex A 
the Hamiltonian is no longer selfadjoint. This means that 
the eigenfunctions are no longer orthogonal. We rather 
have a biorthogonal system which can be normalized as 



(&(A)hMA)> = Ji 



(6) 



where \ipk) and {ipk\ are the right hand and left hand 
eigenvectors of H, respectively. Note that Eq.(|6|) causes 
problems at the exceptional point, since it is exactly at 
this point where two linearly independent eigenfuctions 
no longer exist. This is in contrast to a genuine degen- 
eracy of a selfadjoint operator where a fc-fold degeneracy 
always gives rise to a fc-dimensional eigenspace. At the 
exceptional point, not only the eigenvalues but also the 
eigenfunctions coalesce. As a consequence, the orthog- 
onality conflicts with the normalization, in other words, 
if Eq.(^) is enforced globally, that is also at A = A c , the 
components of the wave function have to blow up. This 
can be made explicit by parametrizing the wave functions 
by the complex angle 9, viz. 
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with 



tan 2 (9(A) 



ffi(A) - E 2 (X) - (ei - eg) - Xjuh - ^ 2 ) cos2</> 
Ei(\) - E 2 (X) + (ei - e 2 ) + A(wi - w 2 )cos2</)' 

(8) 



At A = A c it is Ei = _E 2 and hence tan 2 9 = — 1 which 
implies |cos#| = |sin#| = oo, that is the components of 
the wave functions blow up. (Note that tan 2 9 = in the 
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trivial case <fi — 0). The increase of the components of the 
wave functions while approaching exceptional points has 
been used in similar context as a theoretical signature of 
a phase transition ||, but we do not believe that it has 
observational consequences. 

We now study the behaviour of the wave functions in 
more detail for two contours in the A-plane which start, 
say, at A = and end at large real values of A, but enclose 
an exceptional point between the two contours. For the 
complex angle 9 we choose an expression which is more 
convenient for this purpose, viz. 



tan6»(A) 



X(ll>i — LU2) sin 20 



£i(A) - E 2 (X)+e 1 - e 2 + A(^i -u 2 )cos20 

(9) 



The first path can be taken along the real A-axis. From 
Eq.(^|) we read off the expected result that 9(0) = and 
0(A) — ► (j) for A 3> I (ei — e 2 )/(wi — w 2 )|. In obtaining 
this result use is made of E± — E 2 — 2R — > X(u>i — oj 2 ) 
for A > I (ei — e 2 )/(wi — w 2 )|. For the second path we 
move into the upper A-plane in order to pass above the 
exceptional point before returning down to the real axis 
again. Using again Eq.(^) we now have to observe that 
we crossed into the other sheet which means E\ — E2 — 
— 2R — > — \(u>i — w 2 ). As a consequence we find this 
time tan# = — cot0 = tan(0 + tt/2). Surely we would 
expect the wave functions to interchange just like the 
energies, if an exceptional point is encircled. But our 
finding indicates that one wave function has changed its 
sign. In fact, we obtain along the second path 



ipi->ip2, i>2 



(10) 



Note that this result is equally obtained, if a closed con- 
tour which encloses the exceptional point is described in 
the A-plane. Of interest is the result of a double loop in 
the A-plane. It yields 



4<i 



-Ipl, 4>2 



~lp2 



(11) 



which is, in accordance with the corresponding single 
loop in the energy plane, just Berry's phase retrieved. 
We mention that this implies that the wave functions, if 
parametrized as in Eq.([7), have an algebraic singularity 
that is determined by a fourth root at the exceptional 
point; only a four- fold loop in the A-plane restores com- 
pletely the original situation. 

Next we address the question concerning the physi- 
cal significance of these findings. It boils down to the 
problem of varying a complex interaction parameter in 
the laboratory. Problems of a similar nature have been 
discussed in connection with Berry's phase for dissipa- 
tive systems |l2| ], but not as yet for exceptional points. 
We are guided by the phenomenological description of 
open quantum systems Q and submit as one suggestion 
a strongly absorptive system, where the parameter A in 



Eq.([j]) is traditionally replaced by — iG with real absorp- 
tion parameter G. With the replacement the eigenvalues 
acquire imaginary parts which are related to the inverse 
life times of the states of the open system. The excep- 
tional points appear now in the complex G-plane at 



G r = - 



ei - £2 



exp(±2z</> + i7r/2). 



(12) 



U>1 - U! 2 

If the coupling of the two channels is equal, that is if 
4> = 7r/4, the two exceptional points lie on the real G- 
axis at 



G c 



± 



ei - £2 



If the coupling is nearly equal, they lie just above or be- 
low the real axis depending on the coupling being slightly 
weaker or stronger (0 < ir/4 or <f> > ir/4). Controlling 
the absorption parameter G and the relative coupling en- 
ables one to pass the exceptional point on a path below 
and above, which is the situation described above. A 
setup similar to the one proposed by Berry Q and im- 
plemented in |ljj should make feasible a distinction of 
the two cases. In a two dimensional electromagnetic res- 
onator the relative coupling of two suitable levels, that is 
levels that do coalesce, can be controlled by a judicious 
choice of the geometry of the resonator The global 
absorption can be controlled by radiation losses regulated 
by suitable antennas. In this way the change of the phase 
of one but not the other wave function should be observ- 
able. Note that, for large values of G, one state will be 
much broader than the other || . 

We stress again that the occurrence of exceptional 
points is a generic mathematical feature associated with 
any system that has avoided level crossing. While this 
presentation used a two dimensional illustration, the re- 
sults can be immediately generalised. In fact, if more 
than one exceptional point is encircled, the resulting 
phase change is simply a combination of several two- 
dimensional cases. The findings presented in this paper 
have therefore universal significance. Whether the par- 
ticular phase changes of the respective wave functions 
are detectable, is subject to experimentation. Due to the 
short life time involved for one of the states, scattering 
experiments, where a steady incoming flux can be main- 
tained, are expected to be particularly suitable. Fortu- 
itously, in the spirit of the experiment performed by |l3| ] 
the adiabatic change of the parameters is not essential, 
since it is the movement of the nodes of the wave func- 
tions that are being observed under parameter variation 
flTij . Possible jumps of the phases do therefore not affect 
the result of interest. 
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